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ABSTRACT 

By assuming an aspherical stellar wind with an equatorial disk from a red giant, we 
investigate the production of Type la supernovae (SNe la) via symbiotic channel. 
We estimate that the Galactic birthrate of SNe la via symbiotic channel is between 
1.03 X 10~^ and 2.27 x 10~^ yr~^, the delay time of SNe la has wide range from ~ 
0.07 to 5 Gyr. The results are greatly affected by the outflow velocity and mass-loss 
rate of the equatorial disk. Using our model, we discuss the progenitors of SN 2002ic 
and SN 2006X. 
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1 INTRODUCTION 

Type la supernovae (SNe la) are exploding stars which 
are good cosmological distance indicators and have been 
used to measure the accelerated expansio n of the Universe 
l|Riess et al.lll998l : iPerlmutter el all 1 19991 ). It is widely ac- 
cepted that progenitors of SNe la are mass accreting carbon- 
oxygen (CO) white dwarfs (WDs) and they explode as a 
SN la when the ir masses reach app roximately the Chan- 
drasekhar mass (jNomoto et al.l | l984| ). Two families of pro- 
genitor models have been proposed: the double-degenerate 
model and the single-degenerate model. For the double 
degenerate model, previous works indicated that the ex- 
pected accretion rates may cause the accretion-induced col- 
lapse of the CO WDs and the formation o f neutron stars 
(|Nomoto fc Ibenll 19851 : ISaio fc Nomotolll985l ') . In the single- 
degenerate model the mass donor is a main sequence (MS) 
referred to as the WD-I-MS cha nnel or a red giant (RG) re- 
ferred to as WD-I-RG chan nel ivan den Heuvel et al.lll992l : 
I Li fc van den Heuvel|[l997l ). The latter is also called a sym- 
biotic channel. 

Accord i ng t o iHan fc Podsiadlowskil (|2004l ') and 
iMeng et all (|2008l), the birth rate of SNe la via the WD-(-MS 
channel can only account for u p to 1/3 of 3 — 4 x 10~'^yr~^ 
observ ation ally estimated bv [van den Bergh fc TammannI 
l|l99ll ) and ICappellaro fc Turattd (Il997l'). Other channels 
should be important. Recentlv. IWang et all l|2009l ) showed 
that heUum star donor channel is noteworthy for producing 
SNe la (~ 1.2 x 10~^ yr~^). Howeve r, helium associated 
with SNe la fails to be detected l|Marion et al.l l2003l : 
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iMattila et "allbOOSh . and iKato et al.l (|2008l ') suggested that 
SNe la from helium star donor channel are rare. Therefore, 
symbiotic channel is a possible candidate. By detecting Na 
I abso rption lines with low expansion velocities, iPatat et al.l 
(|2007l ) suggested that the companion of th e progen i tor o f 
the SN 2006X may be an early RG (However, [Chugail (120081 ) 
showed that the absorption line s detected in SN 2006X 
cannot form in the RG wind). IVoss fc Neleman^ (|2008| ) 
studied the pre-supernova archival X-ray images at the 
position of the recent SN 2007on, and they considered that 
its progenitor may be a symbiotic binary. Unfortunately, 
a WD+RG binary system usually undergoes a common 
envelope phase when RG overflows its Roche lobe in 
popular theoretical model. WD-I-R G binaries are unl i kely t o 
become a main cha n nel for S Ne la. Yungelson et. al.] (Il995h . 
LYungclson & Livio (1998), Han & Podsiadlowskil (|2004l ) 
and Lii ct al. (2006) showed that the birthrate of SNe la via 
symbiotic channel are much lower than that from WD-(-MS 
channel. In order to stabilize the mass tran s fer pr ocess and 
avoid the common envelope, iHachisu et al.l (|l999l ) assumed 
a mass-stripping model in which a wind from the WDs 
strips mass from the RG. They obtained a high birth rate 
(~ 0.002 yr"^) of SNe la coming from WD-I-RG channel. 

Most of previous theoretical works assumed that the 
cool giant in symbiotic stars shares the same mass-loss 
rate with fi eld giant and the wind from the cool gian t 
is spherical (|Yungelson et. al.l Il995l : iLii et al.l [20061 . |2008| ). 
However, based on cm and mm/submm radio o bservations 
JSeaauist et al.lll993l: iMikolaiewska et al.|[2003l ') and IRAS 
data ( Kenvon et al.l 19881 ). mass-loss rates for the symbi- 
otic giants are systemat ically higher th a n tho se reported 
for field giants. Recently, IZamanov et all (|2008l ) found that 
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the rotational velocities of the giants in symbiotic stars are 
1.5 — 4 times faster than those of field giants. Using the rela- 
tion betw een rota tional velo c ity an d mass-loss rate found by 
iNieuwenhuiizen fc de Jageij (|l988h . IZamanov et all (|2008| ) 
estimated that the mass-loss rates of the symbiotic giants 
are 3 — 30 times higher th an those of field giants. In addi- 
tion, lO-Brien et aLri|2006h suggested that the bipolarity in 
the 2006 outburst of the recurrent nova RS Oph may result 
from an equatorial enhancement in its cool giant. If the cool 
giant in symbiotic star has high mass-loss rate and an as- 
pherical stellar wind, the contribution of symbiotic channel 
to total SNe la may be significantly enhanced. 

In this work, assuming an aspherical stellar wind with 
equatorial disk around the symbiotic giant, we show an al- 
ternative symbiotic channel to SNe la. In § 2 we present our 
assumptions and describe some details of the modeling al- 
gorithm. In § 3 we show WD-I-RG systems in which SNe la 
are expected. The population synthesis' results are given in 
S 4. Main conclusions are in 5 5. 



2 PROGENITOR MODEL 

For the binar y evolution, we use a rapid binary star evo- 
lution code of iHurlev et al] (|2002l ). When the primary has 
become a CO WD and the secondary just evolves into first 
giant branch (FGB) or asymptotic giant branch (AGB), an 
aspherical wind from the secondary is considered. In other 
phases of bi n ary e volution, we adopt the descriptions of 
iHurlev et al] l|2002l ). 



2.1 Aspherical winds from cool giants in 
symbiotic stars 

In general, the stellar wind from a normal KG is expected 
to be largely spherical due to the spherical stellar surface 
and isotropic radiation. However, the majority (> 80%) of 
observed plan etary nebulae are found to have as pherical 
morphologies (IZuckerman fc Ailed [l986l : iBalick' 1987). This 
property can be explained by two models: (i)the generalized 
wind-blown bubble in which a fast tenuous wind is blown 
into a previously ejected slow wind (see a review bv iFranl] 
l|l999l )l: (ii)the interaction of the slow wind blown by an 
AG B star with a coUimated fast wind blown by its compan- 
ion l|Sokej|2000l ). In two models, the slow wind is aspheri- 
cally distributed and the densest in the equatorial plane. The 
aspherical wind with an equatorial disk may result from the 
stellar rotation (Biorkman & Cassinolli 1993 ) or a simple 
dipole magnetic field (M att et al.ll2000i ). lAsida fc TuchmanI 
l|l995l ) showed that an anisotropic wind from AGB star can 
be caused by combining rotation effects with the existence 
of an infiated atmosphere formed by the stellar pulsation. 
Therefore, the fast rotation and strong magnetic field are 
crucial physical conditions for an equatorial disk. 

An isolated KG usually has a low rotational velocity so 
that its effect can be neglected. Furthermore, using a tra- 
dition al method applied for the solar magnetic field, ISokeij 
l|l992t) estimated the magnetic activity of AGB stars, and 
found that the level of activity expected from single AGB 
stars is too low to explain the aspherical wi nd. However , 
the situations in symbiotic stars are different. ISoked (|2002l ) 
showed that the cool companions in symbiotic systems are 



likely to rotate much faster than isolated RGs due to accre- 
tion, tidal interaction and back-flowing material. Accord- 
ing to measurements of the projected rotational velocities 
of the cool giants 9 symbiotic stars and 28 field giants, 
IZamanov et al.l (|2008l ) found that the rotational velocities 
of the giants in symbiotic stars are 1.5-4 times faster than 
t hose of field giants, which confirmed the results of ISokej 
(200?). Stellar magnetic activity is closely related to the 
rotation. The cool giants with fast rotational velociti es are 
prim e candidates for possessing strong magnetic field (ISokeil 
l2002h . Therefore, the cool giants in symbiotic stars, having 
the fast rotation velocity and strong magnetic field, can re- 
sult in aspherical winds. However, in binary systems, the 
orbital motion, the magnetic field and the gravitational in- 
fiuence of t he companions can complicate the structure o f 
the outfiow (|Sokeij|l994l . Il997l : iFrankowski fc Tvlendall200ll ) . 
A detailed structure of the winds is beyond the scope of this 
work. By assuming several parameters, we construct a pri- 
mary aspherical model to describe the morphology of winds 
from cool giants in symbiotic stars. 

No comprehensive theory of mass loss for RGs exists at 
present. IHurlev et a l. (2000) applied the mass loss laws of 
iReimersI lfl975h and IVassihadis fc WoodI l|l993h to describe 
the mass-loss rates of FGB and AGB stars, respectively. As 
mentioned in Introduction, the mass-loss rates for the sym- 
biotic giants are systematically higher than those reported 
for field giants. We use a free parameters C, to represent the 
enhanced times of mass-loss rates during FGB and AGB 
phases for the symbiotic giants by 

Ml = CMrg (1) 

where Mrg is the mass-loss rate of RG in IHurlev et al.l 
l|2000l ). 

In the present work, the winds from cool giants in sym- 
biotic stars fiow out via two ways: an equatorial disk and a 
spherical wind. The total mass-loss rate Mi, is represented 
by 

Ml =M'^ + M'P'' (2) 

where M"^ and M=p'' give the mass-loss rates via the equa- 
torial disk and the spherical wind, respectively. We assume 
that the ratio of the mass-loss rate in the equatorial disk to 
the total mass-loss rate is represented by a free parameter 
77, that is. 

Ml = M''/ri = ArPV(l " v)- (3) 

In this work we make different numerical simulations for a 
wide range of C, and 77. 

The mass-loss rates Ml, M'^ and M^'^'^ are affected by 
the rotational velocities and the magnetic activities of the 
RGs in symbiotic stars. The rotational velocities and the 
magnetic activities depend on the binary evolutionary his- 
tory (Sokcr 20o3). In the present paper we focus on the 
effects of the aspherical wind with an equatorial disk on the 
formation of SNe la. We neglect the different rotational ve- 
locities and magnetic activities which result from the differ- 
ent binary evolutionary history and result in different A/l, 
M*^ and M^'^^ . This work overestimates the contribution of 
the sy mbiotic stars with long orbital periods to SN la. Re- 
cently, ISokeij (|2008l ) suggested the existence of an extended 
zone above the AGB star where parcels of gas do not reach 
the escape velocity. In general, the binary with an AGB 
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star have a long orbital periods. The low outflow velocity of 
stellar wind greatly increases the efficiency of the compan- 
ion accreting stellar wind. In our work, we do not consider 
the extended zone above the AGB star, which results in an 
underestimate the contribution of the symbiotic stars with 
long orbital periods to SN la. Therefore, it is acceptable for 
a primary model to neglect the binary evolutionary history 
and the extended zone above the AGB star. 



2.2 Wind accreting 

In symbiotic stars WDs accrete a fraction of the stellar winds 
from cool companions. In the present paper a stellar wind 
is composed of a spherical wind and an equatorial disk. 
F or the former, accordin g to the classical accretio n formula 
in lBondi fc Hovld (|l944h . iBoffin fc Jorisser] ()l988l l gave the 
mean accretion rate in binaries by 
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where 1 ^ ^ 2 is a parameter (Following iHurlev et al.l 
l|200(]t ). 5w = I), fw is the wind velocity and 
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where a is the semi-major axis of the ellipse, t;orb is the 
orbital velocity and total mass Mt = Mhot + Afcooi- For the 
equatorial disk, we neglect the diffusion of equatorial disk, 
that is, the disk thickness is constant. The accretion rate of 
WD is relative to the crossing area of its gravitational radius 
in orbital plane, which is approximately given by : 



(6) 



where Rq — 2GMwd/{v^ +'''orb) is the gravitational radius 
of WD accretor. If Vorb >> Vw, the accretion rate given by 
Eq. (|6]) may be higher than M'^. It is necessary that the WD 
does not accrete more mass than that lost by the RG. So 
we enforced the condition ^ 0.9M'^. For «orb << Vv, 
we requested is not lower than Rg/{Tva)M'^ . Eq. ((6} 
overestimates the contribution of the symbiotic stars with 
long orbital periods to SN la because we neglect the diffusion 
of equatorial disk. Total mass accretion rate is 



Ma = Afa + 



(7) 



Based on Eqs. Q and (|6}, accretion rate depends 
strongly on the wind outflow velocity which is not readily 
determined. For the spherical winds, we adopt the prescrip- 



tion in iHurlev et al.l (I2OO2I ) , and = y 2/3w '^a^"""' where 

/3w ~ 1/8. In general, the equatorial disk has a lower outflow 
veloc i ty than the spherical wind ()Biorkman fc Cassinellil 
ll993l : lAsida fc Tuchmanlll995l ), and the typical wind veloc- 
ity of field giant is between ~ 5 and 30 km s~^. Due to the 
existence of an extended zone above the AGB star (|Sokerl 
I2OO8I ). the outflow velocity of the equatorial disk can be 
very low. In this work the outflow velocity of the equatorial 
disk Vw is taken as 2, 5 and 10 km s~^ in different numerical 
simulations. 



2.3 Mass transfer rate of Roche lobe overflow 

When a secondary overflows its Roche lobe, we assume that 
there is no equatorial disk for the secondary. At this time 
the mass transfer via the inner Lagrangian point (Li) can 
be dynamically unstable or stable. If the mass ratio of the 
components {q — Mdonor/Af accretor) at the onset of Roche 
lobe overflow is larger than a certain critical value Qc, the 
mass transfer is dynamically unstable and results in the for- 
mation of a common envelope. The issue of the criterion for 
dynamically unstable Roche lobe overflow gc is still open. 
iHieUming fc WebbinkI (|l987l ) did a detailed stu dy of sta- 
bility of ma ss transfer using polytropic models. iHan et aP 
(|200ll . 12002! ) showed that qc depe nds heavily on the as sumed 
mass-transfer efficiency. Recently, IChen fc Hani |2008') stud- 
ied gc for dynamically stable mass transfer from a giant star 
to a main sequence companion. They found that gc almost 
linearly increases with the amount of the mass and angular 
momentum lost during mass transfer. In this work, for nor- 
mal main sequence stars gc ~ 3.0 while gc = 4.0 when the 



secondary is in Hertzsprung gap (|Hurlev et aTl 20021). Bas e 
on the polytrop ii 
I WebbinkI l| 19881 ) 1 



on the polytrop ic models in hiiellming fc WebbinkI ()l987l ). 

gave gc for red giants by 

1 



gc = 0.362 + ■ 



(8) 



3x (1-Mc/A/)' 

where Mc and M are core mass and m ass of the donor, 
respectively. In our work, we adopt the gc of I WebbinkI (| 19881 ). 

When g < gc, the binary system undergoes a stable 
Roche lobe mass transfer and the mass transfer rate is cal- 
culated by 



AfL 



:3.O-''[ln(i?d/i?Ld)]'A/0yr- 



(9) 



where J?d and iiLd are donor's radi us and Roche lobe ra dius, 
respectively. Details can be seen in iHurlev et al] (|2002l ). 

When g > gc, the binaries in which the giants overflow 
the Roche lobe immediately evolve to the common envelope 
phase, while the binaries for the main sequence stars over- 
flowing the Roche lobes will undergo the dynamically stable 
mass transfer. 



2.4 Evolution of accreting WD 

Instead of calculating the effects of accret ion on to the 
WD e xplicitly, we adopt the prescription of iHachisu et al.l 
(|l999l ) for the growth of the mass of a GO WD by ac- 
cre tion of hydrogen-rich materi a l from its compani on (also 
see lHan fc Podsiadlowskl l|2004h : iMeng et al.l (|2008l )). If the 
mass accretion rate A/a exceeds a critical value, Afcr, the 
accreted hydrogen burns steadily on the surface of the WD 
at the rate of Afcr. The unprocessed matter is assumed to 
be lost from the systems as an optically thick wind at a rate 
of Afwind = Afa - A^cr (jHachisu et al.ll 19961 ). The critical 
mass-accretion rate is given by 

Afcr = 5.3 X 10''^ ^''^ Z'^ (MwD 



X 



O.4)Af0yr~^ 



(10) 



where X is the hydrogen mass fraction and is 0.7 in this 
work. If the mass- accretion rate Afa is less than Afcr but 
higher than Afst ~ ^Ma, it is assumed that there is no mass 
loss and hydrogen-shell burning is steady. If Afa is between 
iAfcr and |Afcr, the accreting WD undergoes very weak 
hydrogen- shell flashes, where we assume that the processed 
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mass can be retained. If Ma is lower than |Mcr, hydrogen- 
shell flashes are so strong that no mass can be accumulated 
by the accreting WD. The growth rate of the mass of the 
helium layer on top of the CO WD can be written as 



where 



1, 
0, 



Ma 
Mc: 
Ma 



> M„, 

^ Ma ^ |M„, 

< |M„. 



(11) 



(12) 



When the mass of the helium layer reaches a certain 
value, helium is possible ignited. If helium-shell flashes oc- 
cur, a part of the envelope mass is assumed to be blown off. 
The mass accum ulation efficiency for he lium-shell flashes, 
?7He, is given by iKato fc HachisiJ (|2004l ). Then, the mass 
growth rate of the CO WD, Mvvd, is 



MwD = rju^HcM^ 



(13) 



When the mass of accreting CO WD reaches 1.378 Mq, 
explodes as a SN la. 



3 WD+RG SYSTEMS IN WHICH SNE lA ARE 
EXPECTED 

In this section, according to the assumptions in the above 
section, we simulate the evolutions of the binary systems 
with a CO WD and a MS. All input parameters are the 
same with those in case 1, that is, the outflow velocity of 
equatorial disk = 5 km s~^, rj = 0.9 and C, = 10. The 
initial masses of WDs are 0.76, 0.8, 0.9, 1.0, 1.1, and 1.2 
Mq, respectively. The initial orbital periods are from 0.8 
to 11000 days with AlogP' = 0.1 days. The initial masses 
of MSs are between 1.2 and 2.0 M© with AM^ = O.OSAf© 
when M^vD = O.76M0 and M^^ = O.8OM0, between 1.15 



and 3.0 Mq with AMj = O.IM© when M^d = O.gM© 
and M^vD = l.OM©, between 1.1 and 3.5 Mq with AM| = 
O.IM© when M^o = I.IM©, between 1.0 and 8 Mq with 
AM;^ = O.25M0 when M^^ = I.2M0, respectively. 

Fig. [1] shows WD-f-RG binary systems in the initial or- 
bital period-secondary masse (log P\ M2) plane. Filled sym- 
bols give the binary systems in which CO WDs explode even- 
tually as SNe la: Filled squares indicate SN la explosions 
during an optically thick wind phase (Ma 5S Ma); filled cir- 
cles SN la explosions during stable hydrogen-shell burning 
(Mcr > Ma ^ ^Mcr); filled triangles SN la explosions dur- 
ing mildly unstable hydrogen-shell burning (|Mcr > A/a ^ 
|Afcr). Crosses show binary systems which undergo common 
envelope evolution during WD-(-RG phases, empty squares 
represent binary systems which experience stable Roche lobe 
overfiows during WD-I-RG phases, empty circles give binary 
systems which are detached WD-I-RG systems. They can 
produce symbiotic phenomena while they can not explode 
as SNe la. The binary systems which are not plotted can 
not evolve to WD-I-RG phases. They either explode as SNe 
la via WD-I-MS channel, or become WD-|-dwarf or helium 
MS systems. 

According to Fig. [T] the progenitors of SNe la via 
symbiotic channel are mainly split into left and right re- 
gions. The progenitors in the left region have short ini- 
tial orbital periods. Table [T] shows an example. The sec- 



Table 1. An example for the progenitor of SN la with a short 
initial orbital period. The initial mass of the CO WD is 0.8 
Af0, the initial mass of the secondary is 1.9 Mq and the initial 
orbital period is 1.99 days. The first column gives the evolu- 
tionary age. Columns 2, and 3 show the masses of CO WD 
and the secondary, respectively. The letters in parentheses of 
column 3 mean the evolutionary phases of the secondary. MS 
represents the main sequence, HG for Hertzsprung gap. The 
forth column shows the orbital period. Column 5 gives the ra- 
tio of the secondary radius to its Roche lobe radius. The last 
column shows the critical mass ratio qc- All input physical 
parameters are same with those in case 1. 

AGE (10«yr) Mwd{Mq) M2(Mq) P (Days) R2/RL Qc 



0.0000 


0.80 


1.90(MS) 1.99 


0.373 




1352.0 


0.80 


1.90(HG) 1.99 


0.813 




1357.6 


0.80 


1.90(HG) 1.99 


1.001 


4 


1363.8 


1.26 


0.70(FGB) 2.21 


2.229 


0.83 


1364.1 


~1.378 


0.55 (FGB) 3.25 


2.003 


0.86 



ondary overfiows Roche lobe during Hertzsprung gap. Due to 
Mdonor/A/wD < Qc, the progeultor experiences stable mass 
transfer. During Roche lobe overflows, the great deal mat- 
ter of the secondary has been transferred to the CO WD 
so that the CO WD becomes more massive while the sec- 
ondary mass decreases. When the secondary evolves to FGB 
phase, Afdonor/MwD has been lower than qc. Therefore the 
progenitor undergoes a stable mass transfer in FGB phase 
till CO WD explodes as SN la. In the whole process, there is 
no aspherical wind with equatorial disk from the secondary. 
The progenitors in the right region have long initial orbital 
periods so that the secondaries can evolve into AGB phase. 
Table [5] shows an example. Due to the high mass-loss rate 
of RG and the high accretion rate of CO WD resulting from 
an aspherical wind with equatorial disk, the ratio of mass 
of RG to that of CO WD is lower than qc when secondary 
overfiows its Roche lobe. Therefore, the binary system un- 
dergoes a stable matter transfer until CO WD explodes as 
SN la. In the progenitors with longer initial orbital periods, 
the secondaries never overflow Roche lobe and the accretion 
of CO WDs mainly depends on the aspherical wind with 

equatorial disk. 

Compared with the previous work ( Hachisu et al.|[l999l : 
iHan fc Podsiadlowskill2004l : lMeng et al.ll2008l ). the progeni- 
tors in the present work have longer initial orbital periods 
(up to ~ 10000 days) and wider ranges of initial masses of 
secondaries (from ~ 1.0 to 8 Mq). 



4 POPULATION SYNTHESIS 

We construct a set of models in which we vary different input 
parameters relevant to SNe la produced by aspherical stellar 
wind. Table [3] lists all cases considered in the present work. 



4.1 Parameters for binary evolution and binary 
population synthesis 

In order to investigate the birth rate of SNe la, we carry 
out binary population synthesis via Monte Carlo simulation 
technique. Binary evolution is affected by some uncertain 
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Figure 1. — Final outcomes of the binary evolution in the initial orbital period-secondary mass (log P', Mj) plane of the CO WD+RG 
binary, where P' is the initial orbital period and is the initial mass of the donor star (for different initial WD masses as indicated 
in each panel). Filled squares indicate SN la explosions during an optically thick wind phase (A/a ^ Mcr), filled circles SN la explosions 
during stable hydrogen-shell burning (Mcr > A-fa ^ ^Mcr), filled triangles SN la explosions during mildly unstable hydrogen-shell 
burning (^Mcr > Afa ^ iAfcr)- Crosses show binary systems which undergo common envelope evolution during WD-I-RG phase, empty 
squares represent binary systems which experience stable Roche lobe overflows during WD-I-RG phase, empty circles give binary systems 
which are detached systems during WD-I-RG phase. They can produce symbiotic phenomena while they can not explode as SNe la. The 
numbers along the right y axis in panel 6 only are for Af^j^ = 1.2Af0. Details are in text. 



Table 2. Similar to Table [T] but for the progenitor of SN la 
with a long initial orbital period. CHeB means core helium 
burning. The initial mass of the CO WD is 0.8 Mq, the initial 
mass of the secondary is 1.7 A/q and the initial orbital period 
is 1256 days. 



AGE (lOVr) A/wd(M0) M2{Mq) P (Days) Rz/Rl 



0.0000 


0.80 


1.70(MS) 


1256 


0.005 




1873.6 


0.80 


1.70(HG) 


1256 


0.010 




1901.1 


0.80 


1.70(FGB) 


1256 


0.014 




1977.1 


0.98 


1.33(CHeB) 


795 


0.064 




2108.6 


0.98 


1.30(AGB) 


816 


0.125 




2112.3 


1.10 


1.05(AGB) 


666 


1.001 


1.039 


2113.5 


~1.378 


0.55(AGB) 


890 


1.715 


4.64 



Table 3. Parameters of the models for SNe la via symbiotic 
channel. The parameter in the second column is the out- 
flow velocity of equatorial disk. The parameter ri in column 
three gives the ratio of the mass-loss rate in equatorial disk 
to the total mass-loss rate [See Eq. Q]. The parameter C, in 
the fourth column is the enhanced times of the mass-loss rate 
during FGB and AGB phases for the secondaries [See Eq. JTJ]. 



Cases 


?Jw (km s ^) 


■n 


c 


Case 1 


5 


0.9 


10 


Case 2 


2 


0.9 


10 


Case 3 


10 


0.9 


10 


Case 4 


5 


0.75 


10 


Case 5 


5 


0.5 


10 


Case 6 


5 


0.25 


10 


Case 7 


5 


0.9 


5 


Case 8 


5 


0.9 


30 



input para meters. In this work, the rapid binary star evolu- 
tion code of lHurlev et al.l l|2002l ) is used. If any input param- 
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et er is not spec i ally m entioned it is taken as default value 
in iHurlev et al.l (|2002l '). The metallicity Z=0.02 is adopted. 
We assume that all binaries have initially circular orbits, and 
we follow the evolution of both components by the rapid bi- 
nary evol ution code, includin R the effect of tides on binary 
evolution (|Hurlev et al.ll2002l '). 

For the population synthesis of binary stars, the main 
input model parameters are: (i) the initial mass function 
(IMF) of the primaries; (ii) the mass-ratio distribution of 
the binaries; (iii) the distribution of orbital separations. 

A simple approximation to the IMF of lMiller fc Scalg 

l|l979l ') is used. The primary mass is gene rated using the 
formula suggested bv lEggleton et aXI l| 19891 ) 

0.19X 

' ^ (1 - X)0-75 + 0.032(1 - X)0-25 ' ^^"^l 

where X is a random variable uniformly distributed in the 
range [0,1], and Mi is the primary mass from O.SA/© to 
8M0. 

For the mass-ratio distribution of binary systems, we 
consider only a constant distribution (jMazeh et al.l[l993 ). 



n(g) = l, 0<g<l, 

where q = M2 /Mi . 

The distribution of separations is given by 

log a = 5X -f 1, 



(15) 



(16) 



where X is a random variable uniformly distributed in the 
range [0,1] and a is in Rq. 

In order to investigate the birthrates of SNe 
la, we assume simply a cons t ant star formation 
rate ove r last 15 GvrdHan et al.l Il995l). or a single 
starburst l|Han fc Podsiadlowskil bool : IZhang erahl l2005h . 
In the case of a constant star formation rate, we assume 
that a binary with i ts primary mo r e massive than . 8 Mq 
is formed annu ally dPhillipsI Il989l : lYungelson et al.l 1 1993 : 
iHan et al.l Il995l ). In the case of a single star burst we 
assume a burst producing 2 x 10^ binary systems in which 
the primaries are more massive than 0.8 Mq, which gives a 
statistical error for our Monte Carlo simulation lower than 
5 per cent for SNe la via symbiotic channel in every case in 
Table 13 



4.2 Birthrate of SNe la via symbiotic channel 

Fig. [2] shows the Galactic birthrates of SNe la via sym- 
biotic channel. They are between 1.03 x 10"'^ (case 2) 
and 2.27 x 10~^yr~^ (case 6). Different outflow veloci- 
ties (ii„ — 2, 5, lOkms"^) of equatorial disk in cases 1, 2 
and 3 have a great effects on the birthrates with a fac- 
tor of ~ 40. Different ratios of the mass-loss rate in the 
equatorial disk to total mass-loss rate in cases 1, 4, 5 
and 6, t^s {rj = 0.9,0.75,0.5,0.25), give uncertainty of the 
birthrates within a factor of ~ 9. Different enhanced times 
of the mass- loss rate for the symbiotic giants in cases 1, 
7 and 8, (^ = 10,5,30), have a weak effect within a 
factor of 1.6. The obs ervationally estimated the Galactic 
birt hrate of SNe la by [van den Bergh fc TammannI (|l99ll ) 



and ICappellaro fc Turattol l| 19971 ') 



is 3 - 4 X lO'^'yr-V In 
cases 3, 5 and 6, the contribution of SNe la via sym- 
biotic channel to total SNe la is negligible. However, in 
case 2 with low outflow velocities {v„ = 2 km s~^) and 




9 9 5 

log (t/ yr) 



Figure 2. — The evolution of birthrates of SNe la for a constant 
star-formation rate. The key to the line-styles representing differ- 
ent cases is given in the upper left corner. The numbers along the 
right y axis are only for case 2. 



high rj {rj = 0.90), the contribution is approximately 1/3, 
which is comparable with that via WD-I-MS channel in 
iHan fc Podsiadlowskil l|2004l l. Our results are greatly af- 
fected by the low outflow velocity of equatorial disk and 
its mass-loss rate. 

Fig. [3] displays the evolution of birthrates of SNe la 
for a single star burst of 2 x lO'^ binary systems. Using 
observations of t he ev olution of SNe la rate with redshift, 
Man nucci et al. I l|2006t ) suggested that SNe la have a wide 
range of delay time from < 0.1 Gyr to > 10 Gyr. The 
delay time is defined as the age at the explosion of the 
SN la progenitor from its birth. In the single degener- 
ate scenario, the delay time is closely related to the sec- 
onda ry lifetime and t h us the initial mass of the secondary 
M'i, l|Chen fc Lil l2007l: iH achisu e t al.l [2008bl). In WD-HMS 
chann el, Li fc van den H cuvcl (1997h. lHan fc Podsiadlowskil 
(j2004 ). and iMeng et al.l (12008 ') showed that is between 
2 and 3.5 Mq, which indicates that the range of the delay 
time is from ~ 0.1 Gyr to 1 Gyr. In order to obtain a w ide 
delay time from ~ 0.1 to 10 Gyr. iHachisu et al.l (|2008bl ') as- 
sumed that optically thick winds from the mass-accreting 
CO WD and mass-stripping from the companion star by 
the WD wind. Assuming an aspherical stellar wind with an 
equatorial disk from cool giants in symbiotic stars, we give a 
very wide delay time range from ~ 0.07 Gyr to 5 Gyr. SNe la 
with shorter delay time than 0.1 Gyr result from those pro- 
genitors with more massive M2 than ~ 3.5M0, SNe la with 
longer delay time than 1 Gyr from those progenitors with 
lower M2 than 2Mq. Because the delay time determined 
greatly by the initial mass of the secondary. Fig. O can be 
easily explained by the distribution of the initial masses of 
the secondaries for SNe la in the next subsection. 
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binary systems. The numbers along the right y axis are only for daries for the progenitors of SNe la (the number in every case is 

case 2. normalized to 1 and the width of the bin is 0.3 Mq). The cases 

are indicated in the right top. 




Figure 4. — The distribution of the initial masses of the CO 
WDs for the progenitors of SNe la (the number in every case is 
normalized to 1 and the width of the bin is 0.03 Mq). The cases 
are indicated in the middle top. 



4.3 Distribution of initial parameters 

Fig. |4] shows the distribution of the initial masses of the 
CO WDs that produce ultimately a SN la according to our 
models. There are obviously two peaks. The left peak is at 
about 0.8 Mq, and results mainly from the binary systems 
with short initial periods like that showed in Table [T] and 
with long initial periods like that showed in Table [5] The 
right peak is at about 1.1 Mq, and mainly results from bi- 
nary systems with long initial periods. 

Fig. Ogives the distribution of the initial masses of the 



secondaries for SNe la. The distribution shows three regions. 
The left region is between ~ 1.2Mq (l.OM© for case 2) and 
~ 2.OM0. These SNe la have delay time longer than about 1 
Gyr and correspond to the right peak of Fig. |3] The middle 
region is from ~ 2.QMq to 3.0Mq, and results from the 
progenitors with long orbital periods. These SNe la have 
delay time between ~ 0.5 and 1.0 Gyr and correspond to the 
middle peak of Fig. [S] The right region is more massive mass 
than ~ S.OMq . These SNe la have very short delay time and 
correspond to the left peak of Fig.|3] Except case 6, there is 
a peak at about G.OA/q (also see Fig. [6J. The peak results 
from the binaries in which primary initial masses (MJ) are 
more massive than ~ S.SM©, secondary initial masses are 
between ~ 5.8 and M\, and initial orbital periods are about 
10000 days. For these binaries, before the primaries overflow 
their Roche lobe during AGB phase, the secondaries can 
accrete some material so that the ratio of primary mass to 
secondary mass is lower than the critical value q^- These 
binaries avoid the common envelope evolution. Then, they 
have wide binary separations so that the secondaries can 
evolve RGs and form aspherical stellar winds. 

Fig. |6] shows the distributions of the progenitors of SNe 
la, in the "initial secondary mass - initial orbital period" 
plane. According to Fig. |S] the distribution of the initial 
orbital periods should have double peaks. The evolution of 
progenitors with short and long orbital periods can be ex - 
plained by Tables [T] and respectively. iMeng et all (|2008l ) 
showed that the distribution of progenitors of SNe la via 
WD+MS channel is between ~ 1 days and 10 days. The or- 
bital distribution via symbiotic channel in this work is much 
longer than their distribution. 

4.4 Progenitor of SN 2002ic 

SN 2002ic was the first SN la for which circu mstellar (CS) 
hydrogen has been detected unambiguously (|Hamuv et al.l 
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Figure 6. — Gray-scale maps of secondary initial masses vs. 
initial orbital period P' distribution for the progenitors of SNe 
la. The gradations of gray-scale correspond to the regions where 
the number density of systems is, respectively, within 1 - 1/2, 
1/2 - 1/4, 1/4 - 1/8, 1/8 - of the maximum of Q,,^f/,,^M, , 
and blank regions do not contain any stars. The cases shown in 
particular panels are indicated in their left-up corners. 



I2OO3I '). The evolutionary o rigin o f SN 2002ic has been inves- 
tigated by iLivio fc Reisl (|2003|)_ on the basis of a common 
envelope evolution model, bv lllan fc Podsiadlowskil (|2006l ) 
on the basis of the delaye d dynamical instability m o del o f 
binary mass transfer, by IWood-Vasev fc Sokoloskil (|2006l ) 
on the basis of a re current nova model with a RG, by 
iHachisu et all (|2008al ') on the efficient mass-stripping from 
the companion star by the WD wind. Spectropolarimetry 
observations suggest that SN 2002ic exploded inside a dense, 
clumpy CS environment, q uite possibly with a disk-like ge- 
ometry (jWang et akl 120041 ). Subaru spectroscopic observa- 
tions also provide evidence for an interaction o f SN ejecta 
with a hydrogen-rich aspherical CS medium l|Deng et al.l 
I2OO4I ). 

According to the above descriptions, the progenitor of 
SN 2002ic should have a large amount of CS material which 
has a disk-like geometry. In our model, the mass of the CS 
material is approximately calculated by 



Mcs « -^77Ml + -^(1- 



7?) Ml - AMwD 



(17) 



where D is a distance from SN la, AMwd is the increasing 
mass of WD for a span of and are the outflow 

velocity of the equatorial disk and spherical stellar wind, re- 
spectively. By assuming that S N la explodes inside a spher- 
ically symmetric CS envelope, IChueai fc YungelsonI ()2004l ) 
estimated the total mass of the CS material for SN 2002ic 
is about 0.4 Mq within a radius of about 7 x lO'^^ cm, and 
for SN 1997cy about 6 Mq within a radius of about 2 x 10^*^ 
cm. IWang et al.l (|2004l ) suggested that there are several so- 
lar masses of material asymmetrically arrayed to distances 
of ~ 3 X 10^^ cm. Therefore, we select the progenitor of SN 
2002ic by two conditions: (i)the mass of CS material calcu- 




2 3 4 

Mass of Circumstellar Materials (Mg) 

Figure 7. — The distribution of total masses of CS material for 
SN 2002ic(the number in every case is normalized to 1 and the 
width of the bin is 0.4 Mq). The total masses of CS material is 
give by Eq. I I17I I in which D is 1 X 10^^ cm. 



lated by Eq. (|17|l being larger than O.4M0 within a radius 
of 10^^ cm; (ii)a CS with a disk-like geometry which can 
interact with the ejecta of SN 2002ic. In our model there 
is always an equatorial disk around the secondary before it 
overflows its Roche lobe. 

In Fig. [71 we give the distribution of total masses of CS 
material for SNe la like 2002ic within a radius of 10^^ cm. 
The majority of the CS material lies in the orbital plane. 
Fig. [8] shows gray-scale maps of initial secondary masses M2 
vs. initial orbital period P distribution for the progenitors of 
SN 2002ic. In our model, the mass of CS material depends 
on the initial secondary mass. For examples, there are two 
peaks in Fig. [7] and two regions in Fig. [H] for case 1. The 
left peak around ~ 0.6A'/q in Fig. [7] originates from the low 
region in Fig.[8l and the right peak around ~ G.QMq in Fig. 
[7] from the up region in Fig. [5] 

SN 2002ic appeared t o be a normal SN la from ~ 5 
to 20 days after explosion l|Wood- Vasevll2002l : iHamuv et al.l 
i2003i) . and brightened to twice the luminosity of a normal 
SN la and sh owed strong Hg em ission around 22 days af- 
ter explosion (|Hamuv et al.|[2003l ). The standard brightness 
during the first 20 days after explosion and the suddenness of 
the brightness around 22 days implied that the original SN 
explosion expanded into a region with little CS material, 
and then encountered a region wi th significant CS mate- 
rial (|Wood-Vasev fc Sokoloskill2006l ). Investigating the high- 
resolution optical spectroscopy at 256 d and HK-band in - 
frared photometry at -1-278 and -1-380 d. lKotak et~aLl (|2004l ) 
suggested that there is a dense and slow-moving (~ 100 km 
s~^) outflow, and a dusty CS material in SN 2002ic. In our 
model, a cavity with little CS material around the accret- 
ing WD results from the WD's accretion, the dense and 
slow-moving outflow originates from the dense equatorial 
disk which have collided with the ejecta from SN explosion, 
and dust in CS material is formed in the dense equatorial 
disk. The radius of the cavity is roughly equal to the Roche 
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Figure 8. — Similar to Fig. |6] but for Gray-scale maps of initial 
secondary masses M2 vs. initial orbital period P distribution for 
the progenitors of SN 2002ic. 

lobe radius of the accreting WD, 7?li- The region out of the 
cavity has large amount of CS material which mainly orig- 
inates from the equatorial disk. Due to orbital movement, 
the CS material out of the cavity has complicated struc- 
ture. However, the most dense region lies on the equatorial 
disk around the secondary. The region where the ejecta from 
SN 2002ic firstly encountered the significant CS material is 
between i?Li and binary separation. Fig. |9] shows the distri- 
butions of J?Li and the binary separation of the progenitors 
of SN 2002ic prior to their explosion as a SN la. The obser- 
vationally estimated region is ~ 2 x 10^^ cm away from the 
explosion. Obviously, the majority of binary separations in 
our work are shorter than 2 x 10^^ cm. A few binary separa- 
tions in case 2 in which the equatorial disk has a low outflow 
velocity are longer than 2 x 10^^ cm. As mentioned in i]2.1l 
we do not consider an extended zone above the AGB star, 
which results in the underestimate for contribution of the 
symbiotic stars with long orbital periods to SN la. 

We calculate the birthrates of SNe la like SN 2002ic in 
the Galaxy. From cases 1 to 8, they are 4.0 x 10~^,3.1 x 

0.0, 1.8 X lO"'* and 



^, respectively. Compared with 3— 4x 10 ^ yr ^ 



10"°,2.0x 10~ 
1.6x10"'' yr" 

which is the birthrates of SNe la observationally estimated 
in the Galaxy, the birthrates of SNe la like SN 2002ic in this 
work are not more than 1% of the total birthrates of SNe 
la. SN 2002ic is very rare event in the Galaxy. 

4.5 Progenitor of SN 2006X 

IPatat et all (|2007l ) detected the CS materia l in SN 2006X 
from variable Na I D (However, see IChugail l|2008l )'). They 
found relatively low expansion velocities (Mean velocity is 
about ~ 50 km s~^), and estimated that the absorbing dust 
is a few 10^^ cm from the SN. A possible interpretation given 
bv lPatat et"al] (|2007l ') is that the high- velocity ejecta of nova 
bursts in the progenitors are, by sweeping up the stellar wind 
of the donor star, slowed down to relatively low expansion 
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Figure 9. — Distributions of the Roche lobe radius of the ac- 
creting WD (i?Li) S'lid tliG binary separation of the progenitors 
of SN 2002ic prior to their explosion as a SN la. The left and 
right panels are for Roche lobe radii of the accreting WDs and 
the binary separations, respectively. 



velociti es. According to the resolved light echo, IWang et al.l 
(|2008bl ') fond that the illuminated dust is ~27 — 170 pc. The 
dust su rrounding SN 2006X has multiple shells (|Wang et al.l 
l2008bl ). and stems from the progenitor. In addition, the dust 
surrounding SN 2006X is quite different from that observed 
in the Galaxy, and has a much smaller grain s ize than that 
of typical interstellar dust (|Wang et al.ll2008al lbh. 

Surprisingly, there is apparent presence of dust in the 
CS en vironment in the 2006 outburst of recu rr ent nova RS 
Oph (lO'Brien et al .' 2006; Bode et all |2007| ). lEvans et all 
(^2007^ and iBarrv et al. (.200^ ) reported that the dust ap- 
pears to be present in the intervals of outbursts and is not 
created during the outburst event. This means that there is 
a dense regi on of the wind fro m the RG so th at the dust can 
be for med (iBode et al.l 2007). According to lO'Brien et al] 
(l2006t ) and lBode et all (|2007l ). the densest parts of the red- 
giant wind lie in the equatorial regions along the plane of the 
binary orbit. The dust in recurrent nova RS Oph originates 
possibly from the dense equatorial disk around the RG. 

Therefore, it is possible that the progenitor of SN 2006X 
is very similar with recurrent nova RS Oph. We assume that 
the candidates for the progenitor of SN 2006X satisfy with 
following two conditions: (i) there is a dense equatorial disk 
around RG, which can produce dust; (ii) the progenitor of 
SN 2006X has undergone a series of weak thermonuclear 
outbursts presented by the filled triangles in Fig. [1] prior 
to explosion as SN Icy. The high-velocity ejecta from every 
thermonuclear outburst blow off the parts of the stellar wind 
including the dust produced by dense equatorial disk. Dur- 
ing this process, the high-velocity ejecta are slowed down, 
and the dust can be irradiated so that the grain size becomes 



^ The thermonuclear outburst in this work has an ejecta without 
an increase of CO WD's mass (strong hydrogen-shell burning), 
or has an increase of CO WD's mass without an ej ecta (mildly 
hydro gen-shell burning). However, according to lYaron et al.l 
||2005|) . there is some material ejected during a mildly hydrogen- 
shell burning. We assume that the progenitor of SN 2006X has 
undergone a mildly hydrogen-shell burning. 



10 Lii et al. 



6.8 

4.8 

2.8 

0.. 
6.8 

4, 

2.8 

^0, 
S,6, 

^ 4.8 

2.8 

0, 
6.8 

4.8 

2.8 

0, 



case 1 

1 

- (1) ^ 

— ' \ ' 


cose 2 

1 

. (2) ' ^' 


case 3 


cose 4 

1 

- r . T . 1 - 


case 5 

■ 

^ 

- (5) 


cose 6 

- (6) - 


case 7 

■ 

- . ^ . , - 


case 8 

\ 



2.5 3.0 3.5 4.0 4.5 2.5 3.0 3.5 4.0 4.5 
log P'(days) 

Figure 10. — Similar to Fig. |6] but for progenitors of SN 2006X. 
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Figure 11. — Similar to Fig. |3] but for progenitors of SN 2006X. 



small. A more detailed model for SN 2006X is in prepara- 
tion. 

According to the above conditions, we select the pos- 
sible progenitors of SN 2006X from our sample which are 
showed in Fig. [10] The birthrates of SNe la like SN 2006X 
in the Galaxy from cases 1 to 8 are 9.3 x 10~®,2.5 x 
10"^1.0 X 10"'', 6.3 X 10"^1.2 X 10"®, ~ 0.0,3.8 x 10"® 
and 7.0 x 10"®yr"\ respectively. Therefore, SNe la like SN 
2006X are very rare events. Fig. [TT] displays the evolution of 
birthrates of SNe la like SN 2006X for a single star burst. 



5 CONCLUSION 

By a toy model in which the cool giant in symbiotic star has 
an aspherical stellar wind with an equatorial disk, we inves- 



tigate the production of SNe la via symbiotic channel. We 
estimate that the Galactic birthrate of SNe la via symbiotic 
channel is between 1.03 x 10"^ (case 2) and 2.27 x 10"^ yr"^ 
(case 6), the delay time of SNe la has wide range from ~ 
0.07 to 5 Gyr. The results are greatly affected by the outflow 
velocity and mass-loss rate of the equatorial disk. The pro- 
genitor of SN 2002ic may be a WD-I-RG in which there is a 
dense equatorial disk around the RG. The CS environment 
of SN 2002ic may mainly originate from the dense equatorial 
disk. In the progenitor of SN 2006X, the CO WD has un- 
dergone a series of mildly unstable hydrogen-shell burning, 
and the dust in the CS environment may originate from a 
dense equatorial disk around the RG and is swept out by 
the ejecta of every thermonuclear outburst. 
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